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Abstract. 

Several types of systems were put forward during the past decades to show that 
there exist isospectral systems which are metrically different. One important class 
consists of Laplace Beltrami operators for pairs of flat tori in M" with n > 4. We 
propose that the spectral ambiguity can be resolved by comparing the nodal sequences 
(the numbers of nodal domains of eigenfunctions, arranged by increasing eigenvalues). 
In the case of isospectral flat tori in four dimensions - where a 4-parameters family 
of isospectral pairs is known- we provide heuristic arguments supported by numerical 
simulations to support the conjecture that the isospectrality is resolved by the nodal 
count. Thus - one can count the shape of a drum (if it is designed as a flat torus in 
four dimensions. . . ). 



1. Introduction 

Since M. Kac posed his famous question: "can one hear the shape of a drum" 
[T], the subject of isospectrality appears in many contexts in the physical and the 
mathematical literature. This question can be cast in a more general way by considering 
a Riemannian manifold (with or without boundaries) and the corresponding Laplace- 
Beltrami operator. (Boundary conditions which maintain the self adjoint nature of the 
operator are assumed when necessary.) Kac's question is paraphrased to ask "can one 
deduce the metric of the surface (or the geometry of the boundary) from the spectrum?" . 
Till today, the answer to this question is not known in sufficient detail. An affirmative 
answer is known to hold for several classes of surfaces and domains |21 El El (see also 
a recent review by S. Zelditch |3]). However, this is not always true. One of the first 
examples to the negative is due to J. Milnor who proposed in 1964 two flat tori in M}^, 
which he proved to be isospectral but not isometric Since then, many other pairs 
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of isospectral yet not isometric systems were found. M. E. Fisher considered a discrete 
version of the Laplacian, and gave a few examples of distinct graphs which share the same 
spectrum [7J. A general method for constructing isospectral, non-isometric manifolds 
has been designed by T. Sunada jT^. Sunada's technique applies also to discrete graphs 
and the late Robert Brooks ^T], gave a few examples of families of non-Sunada 
discrete graphs. Sunada's method was also used by Gordon et. al. jH] and Chapman 
[Hj to construct isospectral domains in M^. Other pairs of isospectral domains in 
were proposed in [ini and discussed further in [T3j. Sunada-like quantum graphs were 
presented in jlj. 

Milnor's original work on isospectral flat tori in M}^ induced several investigators to 
find other examples in spaces of lower dimensions. Kneser [T3] constructed an example in 
dimension n = 12, and proved that there exist no such pair in two dimensions. Wolpert 
[TK] showed that all sets of mutually isospectral but non-isometric fiat tori are finite 
at any dimension. The first examples in dimension n = 4 were found by Schiemann 
[12], and later by Earnest and Nipp [T2j. These results were generalized by Conway and 
Sloane [TH], who constructed sets of isospectral pairs of fiat tori in n = 4, 5, 6, and these 
sets depend continuously on several parameters. 

The existence of such a large variety of isospectral pairs, suggests naturally the 
question - what is the additional information necessary to resolve the isospectrality. 
We would like to propose that this information is stored in the sequences of nodal 
counts, defined as follows. Consider only real eigenf unctions of the Laplace-Beltrami 
operator and assign to each eigenfunction the number of its connected domains where 
the eigenfunction does not change its sign (such a domain is a nodal domain). The 
nodal sequence is obtained by arranging the number of nodal domains by the order of 
increasing eigenvalues. Sturm's oscillation theorem in one dimension, and Courant's 
generalization to higher dimensions express the intimate relation between the nodal 
sequence and the corresponding spectrum. However, the information stored in the nodal 
sequence and in the spectrum is not the same, and here we would like to propose that 
the additional knowledge obtained from the nodal sequence can resolve isospectrality. 
We address in particular isospectral fiat tori in four dimensions of the type mentioned 
above. The simplicity of the geometry, together with the rich variety of pairs, make this 
class of systems very convenient, especially here, when the approach is explored for the 
first time. 

The paper is organized as follows. In section ^ we shall summarize some of the 
properties of fiat tori, their spectra and eigenfunctions. The fact that the spectra are 
highly degenerate requires a special choice of the basis set of wave functions for which the 
nodal domains are to be counted. The quantity which signals the difference between the 
isospectral tori is defined in subsection (jSUl). The arguments which lead us to suggest 
that this quantity resolves isospectrality are explained for the families of fiat tori in 
four dimensions JH]; and it is presented in subsection ()3.2|) . A summary and some 
concluding remarks will be given in the last section. 
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2. Flat Tori 



A flat torus is a Riemannian manifold which is a quotient of M" by a lattice of maximum 
rank: T = R"/v4Z", where A = (g^^V ■ ■ g^"^)- Thus, the lattice AZ" is spanned by 
the g^'"^ The reciprocal lattice will be denoted by g^^'-*, and (g^''-' ■ g*-''^) = 6r,s- The 
Gram matrices for the lattice will be denoted by G = A^A, and its reciprocal will be 
denoted for brevity hj Q = G^^ = {A^^){A^^)^ . In the present work we shall deal with 
dimensions n > 4. We shall assume throughout that the lattice vectors A cannot be 
partitioned to mutually orthogonal subsets. 



2.1. Spectra 

We are interested in the spectrum of the Laplace Beltrami operator A — — X^j=i 
with eigenf unctions which are uniquely defined on T. They can be explicitly written 
down as: 

^q(x) = exp||27rifjg,(gW-x)j (1) 

where q = ■ ■ ■ , (Jn) ^ x G T. The corresponding eigenvalues are 

Eq = (27r)2(q.gq). (2) 
The spectrum of a fiat torus may be degenerate, and we denote the degeneracy by 

^7q(^) = Hq e z" : E = Eq} . (3) 

If the matrix elements of Q are rational, it is convenient to express the energy in units 
of (27r)V/(Q) where 1{Q) is the least common denominator of the elements of Q. In 
these units the energy values are integers, and gqi^E) equals the number of times that 
E can be represented as an integer quadratic form. The integer vectors which satisfy 
(j21) will be called representing vectors in the sequel. Their tips are points on the n 
dimensional ellipsoid and their distribution on the ellipsoid will be discussed in the 
sequel. The study of the spectrum (i.e. those integers that can be represented by a 
given integer quadratic form) and the degeneracies (the number of representations) is 
a subject which was studied at length in number theory. Here, we shall give a brief 
summary of the results which are essential for the present work. The interested reader 
is referred to 1201 for further references and details. 

(a) For integer Q, the eigenvalues are integers E which must satisfy the congruence 

E = m{ mod c{Q) ) ; m e N (4) 

and c{Q) is an integer which depends on Q. This result implies that the spectrum is 
periodic. 

(h) The degeneracy gqiE) for integer Q increases as E^~^. This estimate can be derived 
by a simple heuristic argument. The number of integer grid points in a shell of width 
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6E is proportional to E"2 SE/E^. If Q is integer, E must be an integer, and hence, 
the number of values it can take in the interval of interest is 6E. Thus, 

g AE) ~ g^^^ points ^ ^n_^ 
tl possible values 

From (a) and (h) above, it follows that not all integers appear in the spectrum, and 
the gaps are determined by the gaps which appear in the first c{Q) eigenvalues in the 
spectrum. The distribution of the gap sizes may be very complex, and depend delicately 
on Q. Still, one can define the mean gap size, and this value applies for the entire 
spectrum. 

(c) For irrational matrices of the form Q = Qo + aQi where QoiQi are both integer, 
and a is irrational, the mean degeneracy increases more slowly with E, namely, 
g{E) oc Here, the degeneracy class consists of the grid points which satisfy 
both Eq = q.Qoq and Ei = q.Qiq with E = Eq + aEi, and both Eq and Ei integers. 
Their number is proportional to the volume of the intersection of the corresponding 

^ 2 1 _ 

ellipsoid shells E^~6Eq6Ei/{EqEi)2 . The number of spectral values is now the number 
of integer points in the square of size 6Eq6Ei. Thus, in the irrational case, 

girrat{E) OC ^t'^ . (6) 

Figure 1. shows the dependence of the mean value of g{E), averaged over the 
eigenvalues Ei in an interval of width AE centered at E, 

{9{E)) = ^ Yl 9{Ei) . (7) 

v.\Ei-E\</:^E/2 

M stands for the number of eigenvalues in the interval. The linear dependence of {g{E)) 
at sufficiently large E for rational Q and its independence of energy in the irrational 
case are clearly illustrated. 

(d) For integer Q, n > 4 and in the limit of large E, the representing integer vectors of 
E are uniformly distributed on the ellipsoid: Given a well behaved function /(x) on the 
unit ellipsoid S = {x e M" : x ■ Qx = 1}, 

TifFS 5Z /(^)- / /Wdx=0. (8) 
9iE) ^.^^ Je 

2.2. Isospectral flat tori in four dimensions 

To end this section, we shall describe the 4-parameter family of non isometric 
yet isospectral flat tori which is the system we shall consider in most of the 
examples to be discussed in the sequel. This family which was discovered by J. 
H. Conway and N. J. A. Sloane ^H] includes, as a particular case, the previously 
known example of Schiemann Their spectra are given explicitly by the pairs 

of positive definite matrices, which depend on the 4 positive parameters a, b, c, d: 
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Figure 1. Double logarithmic plot of the spectral averaged degeneracy The 
five data sets correspond to: a.) f3 = 2, h) (3 = 3/2, c) (3 — 5/3, d) f3 = 8/5 and e) 
(3 = (\/5 + 1)/2. a)-d) correspond to integer tori and e) to an irrational. The parameter 
P is explained in the text at the end of section I2T^ The dotted line is a linear function 
of E included for comparison. 



1 

12 



^ 12 



/ 9a + b + c + d 
3a — 3b — c + d 
3a + b — 3c — d 

\ 3a — b + c — 3d 

( 9a + 6 + c + (i 
—3a + 36 — c + d 
a + 6 + 3c — d 
V -3a-6 + c + 3(i 



3a — 3b — c^ d 
a + 95 + c + d 
a — 36 + 3c — (i 
a ^ 3b — c — 3d 



3a^b — 3c — d 
a — 36 + 3c — (i 
a + 6 + 9c + d 
a — 6 — 3c + 3(i 



3a — 6 + c 
a + 36 — c 
a — 6 — 3c 
a + 6 + c - 



-3d \ 
- 3d 
-3d 
9d I 



and (9) 



—3a + 36 — c + (i 
a + 96 + c + d 
a + 36 — 3c — d 
a — 36 — c + 3(i 



—3a + 6 + 3c — (i 
a + 36 — 3c — d 
a + 6 + 9c + d 
a — b^ 3c — 3d 



-3a - 6 + c + 3(i \ 
a — 3b — c + 3d 
a — b^ 3c — 3d 
a + 6 + c + 9(i 



Several properties of these matrices can be derived by straightforward computations: 
i. The spectra of and Q~ are identical, and consist of the values a, 6, c, d. The 
unitary matrices which bring to diagonal form are independent of the parameters. 
Explicitly 



D = diag{a, b, c, 4 = T± g± (T^) 



(10) 
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a. and Q~ commute only when at least three of the parameters a, b, c, d are equal. 
From now on when we shall refer to the four parameters family of isospectral tori, we 
shall exclude this set which represents ellipsoids with cylindrical or spherical symmetry. 
Hi. The unitary matrix U which transform to Q~, Q~ = U^Q^U is independent of 
the parameters a, b, c, d, 



U = (T^ 



V 



1 \ 

1 

-1 

-1 / 



det [/ = 1 . 



(11) 



Consider an integer vector q with q • Q^q = E. f/q is on the ellipsoid generated 
by Q~^, but it is integer only if Yli=i \ is even. In this case, also ^^^^ | (Uq)- \ is even. 
Hence, integer vectors with even sums map to each other under U. The integer vectors 
with odd sums do not have this property. 

Conway and Sloane's paper offers another family of isospectral tori in four 
dimensions. This is a two parameter family which we shall not discuss in detail, although 
the numerical results obtained for this case support the conclusions derived from the 
study of the four parameters family. 

To end this section we would like to describe the numerical simulations which 
accompany the subsequent discussions. We have calculated the first 120 million 
eigenvalues for five pairs of isospectral tori (together with the corresponding nodal 
sequences defined in the next section). Four of the pairs have been chosen rational, 
one irrational. The four parameters that define a pair (see eq. ©) are all of the 
form {a,b,c,d) = {a,a/ (3'^,a/ (3^,a/ (3^). Here a just rescales the spectrum, for the 
rational tori it is taken such that the tori are actually integer. For the irrational 
pair it has been set to a = 1. Note, that we present all results as a function 
of E / [detQY^^ which is invariant under rescalings of the spectrum. For the other 
parameter we have chosen the five values (3i = 2, 3/2, 5/3, 8/5, (-\/5 + l)/2. The 
latter value defines the irrational pair {P^ = (y^ + l)/2 is the golden ratio), the 
other parameters are rational approximants to the golden ratio along the Fibonacci 
sequence. Integer tori are obtained by setting = 2^, 2^ x 3"^, 3 x 5^, 2^° (which results 
in detQ(a„ A) = 2^0,316 ^ 320,316 x 5^2,2^^ x 5'^). 



3. Nodal domains and isospectrality 



In this section we shall define the nodal sequences of the flat tori under consideration, 
and will show how they can resolve isospectrality. 



Resolving isospectral "drums" by counting nodal domains 



7 



To define the nodal domains, we consider the real counterparts of (P 

(n \ / n 

27r5^g.(g('-)-x) ; vI/(-) = sin 27r5^g,(gM-x)) (12) 
r=l / V r=l 

and to avoid double counting, we must exclude — q if q is included. 

A convenient representation can be obtained by the transformation 



so that 



and 



Gx ; X e T , yeTy= M"/Z" (13) 



r,s=l 



vl>W(y)=cos(27r(q-y)) ; vl>(^-)(y) = sin (27r(q ■ y)) . (15) 

From this representation it is evident that as long as the lattice vectors cannot be 
reduced to mutually orthogonal subspaces, the real wave functions cannot be expressed 
in a product form. The number of nodal domains is defined by lifting the wave functions 
to M" and counting their nodal domains in the unit cell. The nodal domains form parallel 
strips separated by n — 1 dimensional nodal hyper-planes. The number of nodal domains 
is 

n 

z>[q] =2^|g,| . (16) 

This result can be proved by induction, and we start by assuming that all the qi are 
positive. (fTBjl is trivially true for n = 1, where the nodal manifold are points. In 2-d, 
the nodal manifold are lines which are perpendicular to the direction {qi, q2). There are 
2qi such lines which intersect the unit interval on the yi axis at the points which are the 
nodal manifolds of the 1-d problem. There are additional 2g2 — 1 lines which intersect 
the interval 1 < ?/2 < 1 on the line yi = 1. In total there are 2(gi + ^2) — 1 nodal lines 
in the unit interval and therefore 2(gi + ^2) nodal domains. The same argument can 
now be repeated for any n. The case when some qi are negative can be taken care of 
by a proper reflection. It is also clear that the number of nodal domains of any linear 
combination of \l/q^'*(y) and \l/q ■*(y) (fT2|) depends on q only. 

The strip structure discussed above is very different from the morphology of 
nodal domains in separable systems, which are formed from intersections of locally 
perpendicular hyper-planes. 

The special feature of the fiat tori to be discussed here is that their spectra are 
degenerate, and the degeneracies are maximal when the tori are integer. Any linear 
combination of the eigenf unctions in the degenerate space is itself an eigenfunction, and 
its nodal structure may depend on the particular combination. To be definite, we must 
chose a unique representation of the wave functions, for which the nodal count can be 
defined in an unambiguous way. For this purpose we require that the eigenfunctions are 
presented in the form of a single term and the nodal domains are parallel strips in the 
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unit cell. This definition singles out a well defined basis for each degeneracy space, and 
the number of nodal domains is given by (fT^ . 

3.1. Counting nodal domains 

In a previous work [221; we have shown that counting nodal domains can reveal 
important information concerning the classical dynamics of the system under study. In 
particular, arranging the spectrum by increasing values of the eigenvalues, we studied 
the distribution of the normalized nodal counts = where C>n is the number of nodal 
domains of the n'th wave function. Courant's theorem guaranties that ^„ < 1. The ^ 
distributions have characteristic and universal form if the problem is separable, which 
is distinct from the distribution which prevails when the classical dynamics is chaotic. 

As long as the degeneracy in the spectrum is accidental, one can ignore the 
ambiguity induced in the ordering of the spectrum by the degeneracy. In the present 
case, the degeneracy increases as a power of E, and therefore a different quantity, which 
is independent of the ordering of the wave functions within a degeneracy subspace, is 
called for. 

A natural modification of the normalized nodal counts for a highly degenerate 
spectrum is obtained by considering the quantity 

1 1 " 

We shall refer to vq[E) as the nodal count associated with the degenerate eigenvalue E. 
It coincides with the standard definition for non degenerate cases. The nodal count is 
the tool by which we propose to resolve isospectrality. It is defined as an average over 
the representing vectors on the surface of the ellipsoid E = f\ - Qq. Because of (jH)), the 
mean value of vqiEi) for large E assumes a well defined, and geometrically appealing 
value. 

3.2. Nodal counts and isospectrality 

In this section we shall summarize the evidence we have in support of the conjecture 
that the sequences of nodal counts of tori which are isospectral but not isometric are 
different. Denoting by and Q~ the corresponding pair of Gramm matrices, and by 
E a spectral point, we study the difference 

5v{E) = vq+{E) - vq^{E) . (18) 

Rather than examining (fT^ for individual eigenvalues, we shall consider its average over 
spectral intervals, and the fluctuations about the mean. 

Using the fact that in the limit of large E the representation vectors are distributed 
homogeneously over the ellipse ((HI), we get, 

{vq{E)) ocE"^ [ ds (5(1 - s ■ gs)z/[s] 



Resolving isospectral "drums" by counting nodal domains 9 

- ^* /. r'^r^. ■ (IS) 

where the factor is due to the hnear dependence of z>[q] on q. We use this expression 
to prove that {5v{E)) = 0. This is trivially true for a = b = c = d, where the two 
ellipsoids degenerate to a sphere. In [Appendix A we show this in general, i.e. the 



smooth nodal counts agree. The numerical simulations reproduce this behavior, and we 
use this fact as a check of the numerical procedure. 

The fluctuations in the difference of the nodal sequences are best demonstrated by 
studying the variance 

{(HE)?) = ^ 

m 

) - UQ-{Em)f ; Em in the spectrum. (20) 

\E-Er„\<A/2 

Here A < E and M = \\{Em : \E - Em\ < A/2}. The variance {{Su{E)Y) was 
computed for the examples of non isometric yet isospectral tori discussed in conjunction 
with figure 1. The results are shown in Figure 2. The variance does not vanish, and its 
E dependence clearly distinguishes between the rational and irrational pairs. It shows 
the variance for the isospectral tori which were introduced at the end of section ()2.2|) . 
The two extreme values {P = 2 and P = (a/5 + l)/2) show the expected dependence on 
E over the entire range of E values. The data for the intermediate sets of parameters 
follow the expected asymptotic behavior only after E is large enough to "distinguish" 
that the coefficients are rational. The same behavior was also observed for the two 
parameter family of isospectral tori in four dimensions mentioned in 

We are unable to prove that the variance ()2()j) differs from zero, and that its E 
dependence follows the behavior expected from the numerical simulations. However, we 
can provide heuristic arguments which explain the systematic trends observed in the 
data. The underlying assumptions are that 

/. The number of nodal domains for individual representing vectors P[q] = 2 ^"^^^ Ig^l 
(fTB|) fluctuate independently about their mean (fTUj) . 

//. The fluctuations in Suq+^E) and Suq-^E) are also independent. (This is supported 
by the observation made in point Hi. of the preceding section, that only a fraction of 
the representation vectors of the two Gram matrices can be transformed to each other 
by a single rotation.) 

The first assumption together with the central limit theorem leads to the conclusion 
that 

{{uq{E) - ME))f) cx ^^-^ . (21) 

The factor E in the numerator comes because the nodal count scales as E^. The 
denominator (qqIE)) is the mean number of integer points on the ellipsoid. This 
behavior is consistent with recent number theoretical estimates of the rate of convergence 
of the fiuctuations of the distribution of representations to the uniform limit j21|- Using 
the second assumption, we conclude that 
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Figure 2. Double logarithmic plot of the variance {6v{E)'^) (1201 1. The parameters 
are the same as in figure 1. The dashed line is a linear function of E included for 
comparison. 



Thus, for rational where {gQ{E)) oc E, the variance is independent of E, whereas 
for irrational the variance is expected to vary hnearly with E. This behavior is 
indeed observed in the numerical simulations. 

Another numerical test which supports the validity of the above analysis, consisted 
in computing the fluctuations of the integrated nodal counts 



where J\f{E) ~ E"^ is the spectral counting function Af{E) = J2e<E 9Qi^)- '^^^ difference 
between integrated nodal counts for the isospectral pair can be written as 




AN{E) = NQ^iE) - NQ-iE) = | > ^ z>[q] - > ^ z>[q] | , (24) 

where V± := {q G Z : q ■ Q±q < E and q ■ Q=pq > E}. Here, we made use of 
the exact cancellations of the contributions of q G Z in the intersection of the two 
ellipsoids q ■ Q+q < E and q • Q-q < E. There are further exact cancellations since 
i>[q] is invariant under reflections qi —qi and due to the symmetries between the 
two ellipsoids. Assuming that uncorrelated contributions to AN{E) stem from a thin 
layer at the surface of the two ellipsoids while bulk contributions cancel we have a sum 
over s ~ |q|^ = E^^'^ independent contributions of order z/[q] ~ E^^"^ that vanishes in 
the mean. For the variance this assumption leads to {Nq(E)'^) ~ x(e)'^ ~ E~^^'^ in 
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Figure 3. Double logarithmic plot of the squared nodal count difference {AN{E)'^) 
(smoothed by taking an average over a small spectral interval AE). The parameters 
are the same as for figure 1. The dotted line shows the theoretical prediction E^^^'^ . 



accordance with our numerical analysis (see figure 3). Note, that this result is valid for 
both, rational and irrational tori. 

4. Summary and conclusions 

The heuristic arguments as well as the numerical evidence collected above, support our 
conjecture that isospectral flat tori can be distinguished by studying the fluctuations of 
their nodal counts. A formal proof is still lacking. At the same time, the work indicates 
a hitherto unnoticed link between isospectrality and nodal counts. We hope that the 
present work will induce more research effort in this direction. 

The nodal count studied here is a particular function of the representation vectors 
whose fluctuations distinguish between the isospectral tori. Are there any other 
functions which are sensitive to these differences? We investigated this question to 
some extent, and were not able to find a simple criterion which the function should 
satisfy in order to resolve the spectral ambiguity. The nodal count was not chosen 
arbitrarily, and it is rewarding that it does have the required property. 

As was mentioned above, there exist several other known families of flat tori beyond 
the ones discussed here at length. The way of reasoning proposed above should apply to 
these cases as well. Exploratory work in this direction shows that the nodal sequences 
of isospectral pairs are different for the cases that were tried, but a systematic study 
gets prohibitively time consuming as the dimension increases. 

The conjecture that nodal counts resolve isospectrality is now being tested for 
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a different class of operators - the Schrodinger operators on quantum graphs |23j . 
Preliminary results support the validity of the conjecture also for this class of operators 
which are quite different from the flat tori discussed here. 
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Appendix A. The mean nodal counts of isospectral tori 

To demonstrate the equality of the smooth parts consider the simpler quantity 

tq= I da;V|^i|-- (A.l) 

with rj some power e.g. 5/2. The goal is to show Tq+ = Tq-. From the diagonalization 

m- 

= {T^YDT^ (A.2) 

make a change of variables 

o = T^-uj. (A.3) 

This preserves the measure as are orthogonal and the denominators are seen to agree 
as D is the same for the two lattices. The numerator becomes 

By explicit calculation we get for the " + " lattice: 

VT2 \uji\ = |oi + 02 + 03 - 304I + |oi - 02 - 3o3 - 04I (A. 5) 

i 

+ |oi + 3o2 - 03 + 04I + |3oi - 02 + 03 + 04I , 
and similarly for the " — " lattice, 

V12 \uji\ = |oi - 02 + 303 - 04I + |oi - 3o2 - 03 + 04I (A. 6) 

i 

+ I - 3oi - 02 + 03 + 04I + |oi + 02 + 03 + 304I . 

Each absolute value for one lattice can be made to correspond to the other under the 
change of Oj to — Oj. For instance the first term of the " + " lattice matches with the 
fourth of the " — " lattice. In the integral for r such sign changes preserve both the 
integration domain, the measure and the denominator. Therefore when splitting the 
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integral into four, corresponding to the terms of ()A.5IA.6|1 . each of these terms can be 
brought to an agreement. 
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